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1  Objectives: 

Development  of  high  order  methods  for  detonation  waves  and  efficient  Wavelet 
ADI  methods  for  reacting  compressible  viscous  Navier  Stokes  equations  and 
direct  numerical  simulations  of  multi-dimensional  premixed  Flame  Propaga¬ 
tions  and  deflagration  to  detonation  transitions  (DDT). 


2  Status  of  effort: 

We  have  studied  adaptive  wavelet  ADI  method  for  the  2-D  reacting  com¬ 
pressible  Navier-Stokes  equation  for  the  simulation  of  premixed  flames,  and 
high  order  methods  for  2-D  detonation  waves.  During  the  three  year  period 
of  the  grant,  we  have  investigated  the  stability  of  two  dimensional  detonation 
waves,  and  the  role  of  using  wavelet  ADI  as  a  direct  solver  or  as  a  precondi¬ 
tioner  to  achieve  an  efficient  time  integrator  for  implicit  time  integration  of 
Navier  Stokes  equations. 


3  Research  Accomplishments/New  findings 

The  main  accomplishment  for  this  grant  is  a  hybrid  high  order  algorithm 
for  the  study  of  2-D  detonation  waves,  and  the  development  of  an  adaptive 
wavelet  ADI  method  as  a  direct  solver  or  as  a  preconditioner  for  implicit 
time  integration  of  Navier  Stokes  equations.  Convergence  acceleration  using 
wavelet  ADI  as  a  preconditioner  is  studied  for  reaction  diffusion  equations. 
Also  parallel  implementation  of  the  wavelet  ADI  preconditioner  is  investi¬ 
gated. 

Cellulr  structure  of  Detonation  Waves 

Partially  supported  in  this  project,  we  completed  the  development  of  high 
order  hybrid  numerical  simulation  of  two  dimensional  detonation  waves.  The 
major  finding  of  this  work  was  that  the  cellular  structure  of  detonation  waves 
depended  very  sensitively  on  the  numerical  dissipation  of  the  algorithms  rep¬ 
resenting  detonation  fronts.  Further  studies  of  the  work  in  three  dimensional 
cases  were  needed  to  understand  the  three  dimensional  effects  of  detonation 
waves. 
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Figure  1:  Two  dimensional  detonation  wave  Five  snapshots  of  temperature, 
mass  fraction  of  reactant 

Wavelet  ADI  Schemes  for  Multidimensional  Reacting  Navier  Stokes  Equations 

We  have  concentrated  on  the  construction  and  implementation  of  the 
wavelet  ADI  methods  based  on  Beam  Warming  approximate  factorization  of 
viscous  Navier  Stokes  equations  for  two  dimensional  premixed  flames.  The 
wavelet  ADI  methods  can  be  used  as  an  efficient  preconditioner  for  implicit 
time  integration. 

The  3-D  reacting  compressible  Navier-Stokes  Equations  with  a  one-step 
reaction  are  written  in  the  conservation-law  form  as  follows 

du  ,  9F{U) ,  aam  ,  dH(u)  dF.iu)  aa.m  dH,{u) 

(1) 

where  vector  of  conserved  vector  is  17  =  (p,  pu,  pu,  ptu,  e,  pY)^  and  F{U),  G{U), 
H{U)  are  the  inviscid  fluxes  and  Fy{U,  ^U),  Gv{U,  S/U),  Hv{U,  \/U)  are  the 
viscous  fluxes  ,  \jU  =  {Ux,Uy,Uz)  and  =  (0, 0, 0,0,0,  (^(F,  T))^  is  the 
chemical  reaction  source  term  and  T  is  the  temperature  and  Y  is  the  mass 
fraction  of  the  reactant. 

The  Beam- Warming  factorization  scheme  for  (1)  without  the  source  term 
is  based  on  a  second  order  two  level  implicit  time  discretization  and  the 
source  term  can  be  treated  separately  by  a  Strang-type  operator  splitting 
procedure  where  a  stiff  ODE  solver  is  used. 

(l  -h  d.te{DlA^  A  DlB^  +  DlC^))  At/"  =  -AtR^  (2) 

where  Dy  and  D]  denote  the  wavelet  derivative  matrices  along  x,  y  and  ^ 
coordinate  direction,  respectively  and  the  solution  increment  At/"  =  t/""''^  — 


3 


C/"  and  A"  =  and  are  the  Jacobian 

matrices  for  the  flux  vectors,  the  residual  =  Dx{F  —  Fy)'^  +  Dy{G  — 
GyY  +  D,{H  -  HyY  and  6  =  0.5. 

The  ADI  factorization  with  extra  stability  enhancement  terms  for  3-D  for 
(2)  will  be  based  on 

{\  +  eMDlA^ -eiV,){\  F  6/^tDlB^ -eiVy){lAe/^tDlC-^ -eiV,)\U^ 

=  -Ati?"  -h  CeVe  (3) 

where  the  extra  14,14’^  second  order  implicit  numerical  viscosities  and 
Ve  is  an  explicit  numerical  viscosity  to  ensure  the  linear  stability  of  the  ADI 
scheme,  which  are  especially  needed  for  3-D  problems. 

3.1  Relevance  to  the  AF’s  mission  and  Applications 

The  understanding  of  detonation  waves  and  the  process  of  the  deflagration 
to  detonation  transitions  (DDT)  have  military  and  civil  applications  in  the 
following  areas  and  efficient  adaptive  wavelet  algorithms  will  enable  the  de¬ 
velopment  of  simulation  tools  for  those  applications: 

•  Pulsed  Detonation  Engine  (PDE)  powered  propulsion  System  for  mis¬ 
sile  launches.  The  understanding  of  DDT  process  and  the  detonation 
waves  will  produce  insights  on  the  initiation  of  detonation  waves  and 
the  cycle  time  of  detonation  engines  to  achieve  appropriate  cycling  fre¬ 
quency  while  maintaining  high  specific  thrust. 

•  The  safe  handling  of  explosives  in  military  and  industrial  environments. 


4  Personnel  Supported 

Wei  Cai,  faculty,  PI  -  total  2.5  summer  months. 

5  Publications 

The  following  reviewed  publications  are  supported  by  this  AFOSR  project: 

[1]  High-order  Hybrid  Numerical  Simulations  of  Two  Dimensional  Deto¬ 
nation  Waves,  AIAA  J.  Vol.  33,  No.  7,  pp  1248-1255,  (July,  1995). 
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[2]  Direct  numerical  calculation  of  neutral  stability  curve  for  one  dimen¬ 
sional  detonations,  (with  W.H.  Oh,  Y.L.  Zhu),  SIAM  J.  Scientific  Computing 
Vol.  17,  No.  4,  pp.  814-829,  July,  1996. 

[3]  Adaptive  Wavelet  Collocation  Method  for  the  Initial  Value  Boundary 
Problem  of  Nonlinear  PDE’s,  (with  J.Z.  Wang),  SIAM  Numerical  Analysis, 
Volume  33,  Number  3,  June  1996. 

[4]  Cai,  W.,  and  Zhang,  W.,  An  Adaptive  Spline  Wavelet  ADI  method 
for  two  dimensional  reaction  diffusion  equations.  Journal  of  Computational 
Physics,  Vol.  139,  1-35,  (1998). 

[5]  A  fast  wavelet  collocation  method  for  high  speed  circuit  simulation 
(with  Dian  Zhou),  IEEE  transaction  on  circuits  and  systems,  1999. 

[6]  An  adaptive  wavelet  method  for  nonlinear  circuit  simulation  (with 
Dian  Zhou  and  W,  Zhang)  IEEE  transaction  on  CAD,  1999. 


6  Interactions/Transitions 

(a)  Conferences 

(b)  Contact  with  DoD  Personnels 

[1]  Dr.  Michael  Nusca  at  ARC  lab  at  Arberdeen  approving  ground,  MD 
on  the  possibility  of  extending  our  research  on  detonation  waves  and  wavelet 
methods  in  the  study  of  reacting  flows  inside  a  gun  barrel. 

[2]  Dr.  Datta  Gaitonde  (937-255-7127)  at  Wright  Patterson  AFB  regard¬ 
ing  our  results  results  on  the  flame  and  vortex  simulations  with  wavelet  ADI 
methods. 

[3]  Dr.  Jerome  T.  Tzeng,  Weapons  Technology  Directorate,  Army  Re¬ 
search  lab  on  detonation  issues  in  gun  barrels.  Tel.  410-306-0843. 

(c)  Transitions 

(d)  Patent 

None 


7  Honors  and  Awards 
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Flame  and  Vortice  Pair  Interaction: 

Adaptive  Wavelet  Methods 
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Figure  2:  Adaptive  Wavelet  ADI  Methods  for  Reacting  Navier-Stokes  Equa¬ 
tions:  Flame- Vortices  Interaction 
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